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Abstract. 

To investigate the formation and the propagation of relativistic shock waves in viscous gluon 
matter we solve the relativistic Riemann problem using a microscopic parton cascade. We 
demonstrate the transition from ideal to viscous shock waves by varying the shear viscosity to 
entropy density ratio r\l s. Furthermore we compare our results with those obtained by solving 
the relativistic causal dissipative fluid equations of Israel and Stewart (IS), in order to show 
the validity of the IS hydrodynamics. Employing the parton cascade we also investigate the 
formation of Mach shocks induced by a high-energy gluon traversing viscous gluon matter. For 
r\l s = 0.08 a Mach cone structure is observed, whereas the signal smears out for r\l s > 0.32. 



1. Introduction 

Jet quenching has been discovered in heavy-ion collisions at BNL's Relativistic Heavy Ion 
Collider (RHIC). In this context, very exciting jet-associated particle correlations [2] have been 
observed, which indicates the formation of shock waves in the form of Mach cones [3j, induced 
by supersonic partons moving through the quark-gluon plasma (QGP). Measuring the Mach 
cone angle could give us the possibility to extract the equation of state of the QGP. 

Shock waves can only develop in a medium which behaves like a fluid. The large elliptic flow 
coefficient V2 measured at RHIC [1] implies that the QGP created could be a nearly perfect fluid 
with a small viscosity. Calculations of viscous hydrodynamics [5] and microscopic transport 
theory [6l|7j have estimated the shear viscosity to the entropy density ratio t]/s to be less than 
0.4. There is still an open question if this upper limit of the 7]/s ratio is small enough to allow 
the formation of Mach shocks. 

In this work we address the question, whether and when relativistic shock waves and Mach 
cones can develop in viscous gluon matter for given rj/s values. For this purpose we consider first 
the relativistic Riemann problem [8], which we solve within the kinetic theory and the Israel- 
Stewart (IS) theory of viscous hydrodynamics for comparisons. Here we employ the microscopic 
parton cascade BAMPS (Boltzmann Approach of MultiParton Scatterings) [9] and a solver of the 



IS equations, vSHASTA (viscous SHArp and Smooth Transport Algorithm) |10] . Particularly, 
we demonstrate agreements between the two approaches for matter with (extreme) small r]/s 
values and also show deviations when the r//s ratio becomes large, which implies the invalidity 
of the IS theory. Second, we consider a traverse of a high-energy gluon through gluon matter 
and investigate the formation of shocks in form of Mach cones. Preliminary results are obtained 
by using BAMPS. 



2. BAMPS and vSHASTA 

BAMPS is a microscopic transport model solving the Boltzmann equation 

p^'d^f{x,p) = C{x,p) (1) 

for on-shell particles with the collision integral C{x,p). The algorithm for collisions is based on 
the stochastic interpretation of the transition rate [9]. In this study, we consider only binary 
gluon scattering processes with an isotropic cross section, which is adjusted locally at each time 
step to keep a constant r//s value pT | [T2 t [T3]. 

Simulations of particle evolution in space and time are performed in a static box, where the 
whole box is divided into spatial cells with a volume V^eii = AxAyAz. Collisions of particles 
in same cells are simulated by Monte Carlo technique according to the individual collision 
probability within a time step of At, 

p _ ^ At 

Attest Vcell 

where a is the total cross section, v^^ei = (pi +^2)^/(2-^1-^2) denotes the relative velocity of 
the two incoming particles with four momenta pi,P2 and -/Vtcst is the testparticle number. The 
testparticle method is introduced to reduce statistical fluctuations and is implemented such that 
the mean free path is left invariant. 

The relationship between the shear viscosity rj and the total cross section a is given by 
r] = 4e/(15-R*'') [12], where R^^ = 71(^^1'^*'^) = '^^{vj-cia) / 2> is the transport collision rate in case 
of isotropic scattering processes. () stands for ensemble average in local rest frame. We obtain 

V = ^eAmfp (3) 

for an ultrarelativistic massless gas. Here the LRF energy density is e = SnT, n is the LRF 
particle density, T the temperature and Amfp = l/(?i-(wreic)) is the particle mean free path. In 
kinetic equilibrium the LRF entropy density is given by s = 4n — n In A, where X = ti/uq is the 
fugacity, the ratio of the LRF number density to the one in thermal equilibrium no = IGT^/tt^. 

vSHASTA solves the IS equations of dissipative hydrodynamics for shear pressure and heat 
flow. In 1+1 dimensions the relaxation equations for heat conductivity and shear stress are 

Dq^ = l(q^^g-q^)-P^^-r^^-P^^, (4) 
''9 

DtT = — {tTnS - Tt) - Inl - In2 " , (5) 



with 



-Tvr 



= lq'[e,+Dln^^ , (6) 

Ig2 = -q'Vzlz (dtVz + Vzd^Vz) , (7) 

Iq3 = I bz (.'^zdtTT + d,Tr) + -r,TT {vj, + -i^dtv^)] (8) 
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(q'^'i) {dtv, + , (10) 



v.dtq' + d,q' - ^9,j , (11) 

where qf^g and 7rjv5 are the Navier-Stokes values [10] and D = u^d^. Tq and T,r are the relaxation 
times, respectively. For vanishing shear pressure and heat conductivity the IS equations reduce 
to the equations of ideal hydrodynamics. 



3. The relativistic Riemann problem 

The relativistic Riemann problem [8j is a well-known shock problem in ideal hydrodynamics. 
Initially matter is separated by a membrane at z = in two regions, z < and z > 0, with 
two different pressures Pq and P4. The velocities on both sides are = ^"4 = 0. The matter is 
assumed to be homogeneous in transverse direction. 

Removing the membrane in the case for a perfect fluid, i.e., rj/s = 0, we observe (the green curve 
in Fig. [U a propagation of a shock wave to the right with a larger velocity than the speed of 
sound and a rarefaction wave to the left exactly with the speed of sound. 
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Figure 1. Spatial profile of pressure (left panel, normalized by the initial value) and collective 
velocity (right panel) at t = 3.2fm/c. 

The BAMPS solutions for various rj/s are depicted in Fig. [TJ In particular, the BAMPS 
result for r}/ s = 0.001 reproduces the ideal solution for a perfect fluid to a very high precision. 
With the increasing ry/s value we see a clear transition from the formation of shock waves in ideal 
fluid to the smearing out in free streaming of particles [llj. The characteristic shock profiles 
disappear for large 77/s values. 

In Fig. [2] we show comparisions between BAMPS and vSHASTA for the shear pressure 
TT = 7r^^/7^ and heat flow q^ = h'y'^[vN^ — N^) at t = 3.2fm/c. vr'"^ is the shear stress tensor, 
A^^ is the particle four-flow, 7^ = (1 — v'^)~^ and h = {e + P)/n is the enthalpy per particle. For 
f]/ s = 0.01 we see a very good agreement between vSHASTA and BAMPS, whereas for rj/s = 0.1 
deviations in the region of the shock front and rare faction fan appear. In the case of ry/s = 0.1, 
the local system at the shock front and rare faction fan is strongly out of equilibrium and, thus, 
the applicability of the IS theory is questionable. The microscopic transport approach does not 
suffer from this drawback. 
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Figure 2. Comparison of BAMPS and vSHASTA for the shear pressure tt and heat flow . 



4. Mach Cones in BAMPS 




Figure 3. Spatial profile of the energy density (left panel) and the velocity field (right panel). 
The r]l s value of the thermal medium is l/47r. The initial velocity of the high-energy gluon is 
in z-direction. The length of the arrays in the right panel is unit, i.e., only the direction of the 
velocity is shown. 



Formation of 3-dimensional shock waves in form of Mach cones is investigated by shooting a 
gluon with energy of 20 GeV into a thermal gluon medium with a temperature of T = 400 MeV. 
The thermal medium is embedded in a static box. 

Fig. [3] shows spatial profiles of the energy density and velocity field at t = 2.5fm/c for a 
medium with r]/s = l/An. The energy, which the gluon probe lost due to interactions with the 
medium, creates a shock wave that propagates in form of a Mach cone. The energy density of 
the region behind the Mach cone is smaller than the initial energy density of the medium. This 
region is called a diffusion wake. Collective behavior of the medium response is also clearly seen 
in the profile of the velocity field. Our results agree qualitatively with those found in \14:\ 115] , 
For higher values of the r]/s ratio the typical Mach cone structure smears out as observed in Fig. 
m The strength of the Mach cone signal and also the lower energy density region behind the 
shock front become weaker because of weaker particle interactions in medium with larger rj/s. 




5. Summary 

We have solved the relativistic Riemann problem using BAMPS and vSHASTA, in order to 
investigate the formation of shock waves in viscous matter. A transition from a characteristic 
shock profile in ideal fluid to a complete smearing out in particle free streaming has been observed 
when increasing the r//s value. For high r]/s values deviations between BAMPS and vSHASTA 
occur, which indicates the break down of the IS hydrodynamics. We have also investigated the 
formation of Mach cones within BAMPS. For -q/s = 1/4tt a Mach cone is clearly visible. For 
larger values of ry/s = 0.32 and 1.0 the dissipative effect is so strong that the typical Mach cone 
structure disappears. 
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